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OBSERVERS EOR COMPRESSIBLE NAVIER-STOKES EQUATION 

AMIT APTE*, DIDIER AUROUXt, AND MYTHILY RAMASWAMYt§ 


Abstract. We consider a multi-dimensional model of a compressible fluid in a bounded domain. 
We want to estimate the density and velocity of the fluid, based on the observations for only velocity. 
We build an observer exploiting the symmetries of the fluid dynamics laws. Our main result is 
that for the linearised system with full observations of the velocity field, we can find an observer 
which converges to the true state of the system at any desired convergence rate for finitely many 
but arbitrarily large number of Fourier modes. Our one-dimensional numerical results corroborate 
the results for the linearised, fully observed system, and also show similar convergence for the full 
nonlinear system and also for the case when the velocity field is observed only over a subdomain. 
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1. Introduction. Data assimilation is the problem of estimating the state of 
a dynamical system described by an evolution equation, typically partial differential 
equations (PDE), using observations, often noisy and partial, of that system. This 
has been widely studied in the geophysical context, e.g. meteorology, oceanography, 
fluid flows, etc., [5, 21, 25]. One of the approaches to this state estimation problem, 
and the one that we study in this paper, is the construction of appropriate observers 
[32, Chap. 7]. 

Observers are essentially a modification of the original evolution equation for 
the system, to incorporate the observations in a feedback term, with the aim that 
the solution of the observer converges to the solution of the original system being 
observed. Observers for hnite dimensional systems have been well studied in the 
literature, see, for example [32, 37]. But in many applications such as earth sciences 
or engineering, the systems are modeled using PDE that are highly nonlinear and, 
in many instances, chaotic. In such cases of infinite dimensional systems governed 
by PDE, there are only a few examples available in the literature, mostly for linear 
systems and very few for nonlinear systems [23, 31, 36, 12, 16, 17, 20, 35, 26, 6, 30]. 

Some of the commonly used observers are Kalman Filters or Luenberger observers 
but the main drawback of these is that they often may break intrinsic properties of the 
model, e.g. symmetries and/or physical constraints such as balances in geophysical 
models, see for example [15] and references therein. For nonlinear system possessing 
certain symmetries, it is natural to seek a correction term which also preserves those 
symmetries. Such an invariance may make the correction term nonlocal but it may 
have other desirable properties, for example, independence from the change of coor¬ 
dinates. Recently, there have been attempts to construct, for a variety of systems, 
observers based on considerations of symmetry [9, 10, 3]. The observers we construct 
are motivated by these recent works, as we will see in section 2. Another motivation 
for observer design is the computational cost: a full Kalman filter is usually too ex- 
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pensive for real applications, due to the size of the gain matrices, while an observer 
might be much more affordable, without degrading the identihcation process. 

The main aim of this work is to develop an observer for a class of PDE inspired 
by earth system applications, all of which use some approximations of Navier-Stokes 
equations for fluid flow. In particular, we propose an observer for the nonlinear PDE 
[equations (2.1)] describing the evolution of a compressible, adiabatic fluid whose 
velocity field is observed either fully or partially. This is justified by the fact that 
velocity observations for fluids have become more available in the last years. We 
can cite floats (e.g. Argo) for oceanography, balloons for meteorology, but also all 
observations generated by optical flow or particle imaging velocimetry from images 
(including satellite images). But as we will see, the other case in which the density 
is observed and not the velocity is very similar. Apart from a purely theoretical 
motivation, the observations of sea surface density are now available [13, and references 
therein] and thus the study of observers with density observations is relevant. 

Our main theoretical result, theorem 3.6, supported by extensive numerical results 
of the one dimensional system (section 5), is that in the case of complete observa¬ 
tions of the velocity, an appropriate choice of parameters leads to convergence of the 
observer to the true solution at any prescribed rate for arbitrarily large but finite 
number of Eourier modes. 

There has been some work on developing observers for Navier-Stokes based sys¬ 
tems. The papers [35, 20, 28] work with finite dimensional approximations of the PDE 
involved whereas we work with the full PDE itself. A completely different approach 
based on developing observers using appropriate continuous time limits of discrete 
time 3D-var or Kalman filter is developed in a series of papers [22, 24, 8] and also 
in [7, 2, 1, 27, 18, 4]. All these papers are in the context of incompressible flows. 
Here we deal with the compressible model, a coupled system for density and velocity. 
The linearized system consists of a hyperbolic and parabolic PDE. This coupling of 
mixed types poses some difficulties in tackling the system theoretically, in particular 
the full nonlinear system, unlike in the case of incompressible Navier-Stokes system. 
We refer to [14] and [11], where these difficulties have been overcome for the nonlinear 
system, with highly involved techniques. In the case of the linearized system around 
constant steady states, we manage using Fourier series. One of the main contributions 
of this paper is that we also derive the decay rates for the convergence of the observers 
and indeed find observers that can decay arbitrarily fast. Since Navier-Stokes based 
PDE are commonly used in practical data assimilation problems in earth sciences, 
our work has the potential to be directly relevant to these applications, as we discuss 
in section 6. 

The outline of the paper is as follows. The multi-dimensional model and the 
observer are both introduced in the next section 2. In that section, we also state 
the linearized version of this problem. We analyze the convergence of this proposed 
observer for the linearized PDE in section 3 and prove the main result theorem 3.6. 
We also briefly discuss the difficulties that arise in the theoretical analysis of the 
nonlinear equation or of the cases with either partial observations or unknown forcing 
for the linearized equations. In section 5, we present one-dimensional numerical results 
to substantiate the linear theory. We also present the numerical results showing the 
efficacy of the observer in estimating the true solution for the partially observed linear 
system, the observer in the case of unknown forcing term, and the fully nonlinear case. 
The last section 6 discusses some future directions of research. 

The authors would like to thank both the anonymous referees for useful remarks 
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which led to this improved revised version. 

2. Compressible Navier-Stokes equations and observers. In this paper, 
we study a model of compressible fluid in a bounded domain. The density p{t, x) 
and velocity u{t, x) of the fluid form the state vector of this model and they obey the 
following compressible Navier-Stokes system in n-dimensions: 

(2.1) pt+y ■ {pu) = 0 , p[ut + {u- V)w] = -Vp{p) + pAu -I- (A -I- /r)V(V • u), 

where A and p are the Lame parameters, satisfying the standard assumptions of /r > 0 
and A -I- 2/r/3 > 0 [34, 33]. The pressure p is given by the adiabatic equation of state, 
p{p) = p'^ with 7 , the adiabatic exponent, taken to be 1.4 for the numerical results 
discussed in section 5. We consider solutions over a finite time interval [0,T] for the 
space domain [0,1]” C K" with periodic boundary conditions: 

(2.2) pit, x) = pit, X + ek), V it, x) € [0, T] X [0,1]" , V 1 < fc < n , 

and similar periodic conditions for uit,x) and for all derivatives of p and u, where 
Cfc = (0,..., 1,..., 0) with the 1 in the fc-th place. 

We assume that the initial conditions p(0,x) = piix) and M(0,a:) = uiix) are 
unknown, but we have some information on the solution u(t, x) of equations (2.1). The 
goal is to build an observer (p, u) for this system, in such a way that both the density 
and the velocity of the observer (p, u) converges towards the solution (p, u). We will 
assume that we have observations of the velocity denoted by u(t, x) for t G [0, T] and 
X G n C [0, Ij" where 11 is a subset of [0, Ij". For most of the theoretical study we will 
consider 17 = [0,1]" while for the numerical study of the one-dimensional problem, we 
will consider 11 = [0, L] with 0 < T < 1. 

2.1. Observers. We introduce the observer ip,u), based on the symmetries for 
the system, satisfying the following set of equations: 

.2 3 n Pt+^ ■ ipu) = Fpip,u,u), 

p [ut + i'll ■ V)m] = -Vp(p) -I- pAu -I- (A -I- p)V(V • u) + Fuip,'u,u). 

As there are no observations of the density, we first assume that the feedback terms 
do not depend on p. Additionally, since these terms should be equal to 0 when u and 
M coincide, it is reasonable to consider the following classes of functions: 

(2-4) Fpip,u,u) = ipp* Dpiu-'ll), Fuip,u,u) = ipu* F)uiu-u), 

where Dp and Du are differential (or integral) operators (e.g. V, dt, ■ ■ ■), and ipp and 
ipu are convolution kernels. We assume that these kernels are time independent and 
hence only functions of x, and we assume that they are isotropic. With this choice, the 
observer preserves the symmetries of the system, since the correction terms are based 
on a convolution product with an isotropic kernel, making it invariant by rotation or 
translation. Moreover, observations may be noisy, so the involvement of a convolution 
kernel leads to an observer that is inherently robust to noise. 

Note that for the case of observations over partial domain, i.e., when fl ^ [0,1]", 
the feedback terms are present only in fl so that we essentially write the feedback 
term as: 


(2.5) Fpip,'u,u) = (fip* Dpliu- 'u)1q] , Fuip,'u,u) = (fu* Du[iu-' u.)1q] . 
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We would like to examine the convergence of the observer solution (p, u) to the 
solution of the observed system (p, u) asymptotically in time. In particular, we are 
interested in the rate of convergence of ||p — p|| and of ||{t — u|| towards zero where we 
will only consider the norm in this paper. 

Probably, the simplest observer is defined by Luenberger observer (or asymptotic 
observer): as only u is observed, then the feedback term is added only in the velocity 
equation {Fp = 0), and the feedback term is simply Fu{p,u,u) = — u), where 

fcu > 0 is a constant. We will study more general observers, with the aim of correcting 
also the density equation with the velocity, and in the process increasing the rate of 
convergence of the observer towards the true solution. 

In order to study the theoretical behavior of the observer, we hrst linearize the 
system around a steady state of the nonlinear system (2.1) and study the convergence 
of the observer for the linearized system. We postpone the presentation of numerical 
results for the nonlinear case to section 5.5. 

2.2. Linearization around an equilibrium state and linear observers. 

We now consider an equilibrium state consisting of constant density and velocity 
(po,ito), and we linearize equations (2.1) around this state: 


pt + {uo ■ V)p + poV • M = 0 , 

po [ut + (uq ■ V)it] = pAu + (A + p)V(V • u) - 7 Pg“Vp, 


where only terms linear in (p, u) appear. Of course linearizations around non-constant 
density and velocity will contain additional terms such as (u ■ V)rto and we will not 
consider this case in this paper. 

Remark 2.1. Note that in dimensions n > 2, stationary solutions of equa¬ 
tions (2.1) could be spatially varying, hut in one-dimensional case n = 1, we have the 
following result: 

Proposition 2.2. The only stationary solutions of ID compressible Navier- 
Stokes equations (2.1) are constants: p{t,x) = po, u(t,x) =uo. 

Proof. Let us write the stationary solution as p{t,x) = po{x) and u{t,x) = uq(x). 
Then the first of equations (2.1) in one dimension gives poix) = C/uq{x) for a con¬ 
stant C determined by the initial conditions. The second of equations (2.1) now reads: 

(^C^Uq ^ -\- Cu(f)x = r'QUQxx ^ '^Ox F Cuqx — 


Let us multiply this equation by uq and integrate over the space domain [0,1]; 


-7C">' 


Uq "^Uqx dx + C 


UqUqx dx = V UqxxUq dx. 


The two first integrals are both equal to zero, by considering an integration by part, 
as the boundary conditions are periodic. Then 


0 = -n {uoxf dx, 


meaning that uqx = 0. □ 

We can rewrite the solution along the characteristics of the equation. The trans¬ 
port coefficient for both density and velocity is uq, hence noting that the derivatives 
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of p(f, X + U[)t) are simply Vp(t, x + u^t) and pt(f, x + u^t) + (ug ■ V)p(t, x + Ugi), and 
similarly for derivatives of u(f, x + Ugt), the equations (2.6) at point (t, x + Ugt) are 


(2.7) 


Pt + PoV • M = 0 , 

= pArt + (A + p)V(V • u) - 7 Pq“ Vp, 


which is nothing else than equations (2.6) with ug = 0. So we can now assume that 
ug = 0 without any loss of generality. Thus we will work with this linearized Navier- 
Stokes system with initial conditions p(0, x) = pi(x) and m(0, x) = uj(x), and periodic 
boundary conditions as in equation (2.2). 

As stated above, we assume that the initial conditions pi(x) and u/(x) are 
unknown, but we have observations of the solution u(t, x) of the linearized equa¬ 
tions (2.7). Again, the goal is to build an observer (p, m) for this system, in such a 
way that the observer (p,u) converges towards the solution (p, m). We will use the 
same observer as introduced in equations (2.3)-(2.5), except that the left hand side is 
now linear, just as in equations (2.7): 


Pt+Po'^ ■u = 'fp*Dp[{u-u)lQ] , 

pgUt = pAu -I- (A -I- p)V(V • u) - 7 Pq“ Vp + ipu* Du [(u - M)la] , 

with periodic boundary conditions, and initial conditions different from the true initial 
conditions: p(0,a::) = pi{x) ^ Pi[x) and u{0,x) = ui[x) ^ uj{x). 

We will now present the main theoretical result about this observer in the linear 
case, stating that we can choose the feedback terms in order to guarantee any specified 
rate of convergence for arbitrarily large but finitely many number of Fourier modes. 
Since we are dealing with linear PDF, we will use Fourier series representation as our 
main tool. 

3. Theoretical study of an observer for n-dimensional linear Navier- 
Stokes system. For the case of velocity observations over the full domain (f2 = 
[0,1]"), subtracting the observer equations (2.8) from the state equations (2.7), we 
get the following equations for the errors r = p — p and v = u — u: 


n + PoV • V = -ifp * DpV , 

pgVt = pAv + {X + p)V(V • v) - 7 p[]“Vr - (p„ * DuV , 
which are exactly identical to equations (2.8), with {r,v) replacing {p,u) and with 

M = 0. 

3.1. Damped wave equation formulation. We now eliminate the density in 
the velocity equation, in order to get an equation for the velocity alone. As we will 
see, this will allow us to more easily define the observers, and more particularly the 
differential operators Dp and Du- Starting from equations (3.1), we take the space 
derivative of the density equation, and the time derivative of the velocity equations: 

Vrj -I- pgV(V ■ v) = -ipp * V{Dpv), 

Povtt = pAvt -I- (A -f p)V(V • vt) - 7p[]“Vr4 - pu * DuVt - 

We replace Vr^ in the second equation by its expression given by the first equation, 
and we obtain: 


(3.2) pgVtt= pAvt + {\ + p)y(y ■Vt) + -iPgV{V ■v)+^pI ^pp*V{Dpv)-pu*DuVt. 


Equation (3.2) is a damped wave equation, with two forcing terms coming from 
the observer feedbacks. In this case, the goal is to make the difference v = u — u 
between the observer u and the observed velocity u converge towards 0. We want to 
choose Dp and Du which respect the symmetries of the original system, but of course 
the choice is not unique. Since the observations may be noisy, we also look for the 
operators with minimum number of derivatives. With this in mind, we now assume 
that the following differential operators are used: 


(3.3) 


Dpif) = PoV • /, Duif) = pof, 


which means that, up to constants, we want the velocity equation to be controlled 
by the velocity, and the density equation by the divergence of the velocity. In the 
damped wave equation (3.2), 

Using this choice, equation (3.2) rewrites as follows: 

(3.4) povtt = pAvt + (A + m) V(V • ut) + 'yp^ViV ■ v) + ')pl<Pp * V(V • u) - poPu * vt , 

from where we can see that the last term (with ipu) controls the first two terms on 
the right while the fourth term (with ipp) controls the third term. 

Using equations (3.3), the observer error equations (3.1) now read: 


rt + PoV • V = -poPp * V • u, 

povt = pAv + (A + /i)V(V • v) - 7Pq'^ Vr - poPu * v. 


(3.5) 


We will refer to this system as the full state observer, as it corrects both velocity and 
density equations. 

3.2. Fourier transform. We first note that equation (3.5) defines the solution 
up to a constant. This, if the mean value m of r(t,x) is not equal to 0, then the 
solution can at best converge towards the constant m instead of 0. Hence we assume 
the following. 

Assumption 3.1. We suppose that the mean value ofr{t,x) is equal to 0, i.e. 
the mean values of the initial conditions of p and p are the same. 

Of course, this is a strong assumption as the initial state could have bias. But 
it is quite often assumed, at least in the data assimilation community, that bias on 
observations could be removed, for instance with reanalysis methods, and that noises 
are zero-mean. Moreover, the linear system (2.6) defines the velocity up to a constant, 
and the observer that we propose is consistant with this. We also note that we will 
be able to control the mean value of the velocity, as shown below. 

As we are on the domain [0,1]" with periodic boundary conditions (see equa¬ 
tion (2.2)), we can consider the following Fourier decomposition of the velocity and 
density: 


(3.6) 




k 


k 


where ak{t) G and bk{t) G M are the time dependent Fourier coefficients, 
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Note that we assume 69 = 0i from assumption 3.1. We denote by (ppk and the 
Fourier coefficients of the (time independent) functions <fp{x) and (fiuix) respectively. 

Substituting (3.6) in (3.5), we obtain the following equations for mode k: for all 
1 < i < n, 

(3.7) 

/ " \ 

Poaikit) = -{p4:Tr'^\kf+PQ<Puk)aik{t)-{X+fi)ATT^ '^ajk{t)kj ^i2'Khk{t)ki, 


where \k\'^ = 
(3.8) 


Kit) = -poi27r(l + ippk)'^ajk{t)kj. 

1=1 


Introducing ykit) = 


( bkit) \ 
aifc(t) 


€ the system (3.7)-(3.8) rewrites as 

\ ^nkit) J 

(3-9) y'kit) = Mkykit), 

where Mk is the following (n + 1 ) x (n + 1 ) complex matrix: 


(3.10) 


Mk = 


0 


04(1 + (Ppk)k'^ 


-Csk -(ci|fcp + iy 9 ufc)/„ - C 2 fc 0 A: 
where /„ is the n x n identity matrix, 0 is the scalar zero, and 


(3.11) 


/idTT^ {X + y) 4 TT^ 2.0 O 

Cl =-, C2 =-, C3 = 7P(j i2tt, C4 = -po*27r . 


Po 


Po 


Remark 3.2. Note that the Fourier analysis can also be done on the damped 
wave equation (3.2). Substituting (3.6) in (3.2), we obtain the following equation for 
mode k: for all 1 < i < n, 


Poa”kit) = -ip4:Tr^\k\^ + potpuk)aikit)-{X + n)4Tr^ | ^a'-fc(t)fcj ] h 

vl=i 


(3.12) 


A 


Let Zk{t) = 


( aik{t) \ 

^nk (t) 

Kkit) 


-7/Oo(l + 12pfe)47r^ I '^ajkit)kj h. 

vl=i j 


then equation (3.12) rewrites as Kit) = NkZkit), 


\ a'nkit) J 

where Nk is the following 2n x 2n matrix: 


(3.13) 


Nk = 


On 

-05(1 + (Ppk)k 0 k -{ci\k\‘^ + ^uk)In - C 2 k®k 








where 0„ is the n x n zero matrix, and where we recall de values of Ci and C 2 and 
introduce C 5 ; 


(3.14) 


Cl = -, C2 

PO 


(A + /i)47r^ 
Po 


C5 = C3 X C4 = 7Po ^471^ , 


3.3. Spectral analysis. We now look at the eigenvalues and eigenvectors of 
matrix For fc ^ 0, Let {kf ,1 < i < n — l}bea basis of the set of vectors 
orthogonal to k in R". 

Proposition 3.3. 

• For k ^ 0, the eigenvalues of matrix are: 

(3.15) 

^dk — (ci|A)| “h ^uk')^ 

- ((ci + C 2 )|fcp + ^Puk) ± \J ((ci + C 2 )|fcp + Pukf - 4c5(1 + Ppk)\k\'^ 
A±fc =--- 


with 

(3.16) 


/r47r^ 

Cl = -, C2 

Po 


(A + /i)47r^ 
Po 


C5 = 7Po ^471^ . 


Xdk has a multiplicity n—1, and the associated eigenvectors are the following 
n — 1 vectors in R"+^ .• 

A+fc and X-k both have a multiplicity 1 , and the associated eigenvectors are 

) and ( • 

• For k = 0, the eigenvalues of Mq are: 

(3.17) Acio = -^+0 = 0- 



Ado 'kiow has multiplicity n, with the associated eigenvectors: 
i < n} bein 
eigenvector 


0 

Ci 


, {ei,l< 


i < n} being a basis o/R", and A+o has multiplicity 1 , with the associated 

1 


The proof is straightforward and simply consists in multiplying these eigenvec¬ 
tors by Mk- Note that in the case k = 0, the eigenvalue A+o cannot be controlled, 
and the associated eigenvector corresponds to constant density solutions. So, from 
assumption 3.1, the projection of r(x,t),v(x,t) along this eigenvector will be zero. 

For fc 0, the discriminant of the quadratic equation for the eigenvalues X±k is 


(3.18) Afc = ((ci -f C 2 )|fcp -f Puk)'^ - 4 c 5(1 -f (ppk)\k\'^ 


and we use the notation y^= i^/—Ak if the discriminant is negative. 

The eigenvectors associated to Xdk correspond to non constant divergence free 
velocity solutions, for which the equation simply consists in a diffusion equation (plus 
the feedback term), the diffusion coefficient being Xdk- Also note that the first com¬ 
ponent of the two eigenvectors associated to X±k is complex, since C 4 = —po*27r. 

We will now relate the eigenvalues of the matrix Mk to the decay rate of the 
solutions yk{t) of equation (3.9). 
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Definition 3.4. dk > 0 is a decay rate of yk(t) if there exists a constant c > 0 
such that |lyfc(<)|| < for any t > 0. 

Of course, we are usually interested in the largest decay rate. In our particular 
situation, the eigenvalues given by equations (3.15) are all negative if ipuk > 0 and 
‘Ppk >0. In such a case, it is straightforward to see that the largest decay rate is 
given by: 

(3.19) dk = — max {5i(A)} 

AeSp(Mfc) 

because there exists c > 0 such that ||?/fe(t)|| < Thus, the largest decay 

rate can be controlled by ipuk and ippk as we will see in next section. 

Remark 3.5. For k ^ Q, the eigenvalues of matrix Nk are the following: 

(3.20) Aofe = 0, Xdk, A±fc, 


where Xdk and A±fe are given by (3.15). 

The new eigenvalue Aofc has a multiplicity n — 1, and Xdk still has a multiplicity 
n — 1. The n — 1 associated eigenvectors are: ( J , with X = Xok or Xdk- 

Note that the solutions corresponding to Xok ore solutions of the damped wave 
equation (3.2), but not solutions of the original system (3.1). 

These eigenvectors correspond to divergence free solutions, either constant in time 
(for X = Xok) or not (for X = Xdk)- 

Finally, A+fc and X-k both have a multiplicity 1, and the eigenvectors are respec- 

AhIO (aIa)' 

Note that the eigenvalue Xok is artificial, and only appears because we took the 
time derivative of the velocity equation in order to eliminate the density. This is the 
only main difference between Mk and Nk spectral analysis. 


3.4. Main result. We now give the main result in this framework: 

Theorem 3.6. We assume that assumption 3.1 holds. For any d > 0, for any 
K > 0, one can find ipp{x) and (pu(x) such that the maximal decay rate of the errors 
r{t,x) and v(t,x), solutions of (3.5), towards 0 is at least d for any Fourier mode 
k such that |fc| < K. The following values can be chosen, with ci, C 2 and C 5 given 
by (3.16); 


(3.21) — d , TpO — 0 ? 

(3.22) = max{0; d-ci|fcp; 2d - (ci + C 2 )|fcn} , 0 < \k\ < K, 

AOOOA i ^ {{ci + C2)\k\^ + (fukf 

(3.23) yjpfe = max0; ^- - -^|’ 0<|fc|<R:, 

(3.24) (^„fc=0, ippk = 0, \k\>K. 


In other words, for any specified decay rate, we can find convolution kernels ipp 
and ipu such that the observer u converges towards u at this specified rate up to 
any Fourier mode. Indeed, we can choose appropriate Fourier coefficients of these 
kernels for any mode k £ IN, but as we will see, their expression does not ensure the 
convergence of the Fourier series and we need to truncate the series. 

Proof. Let d > 0, dsT > 0, and let k a Fourier mode such that 0 < |/c| < K. 
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The case fc = 0 is trivial: the only eigenvalue is leading to the 

appropriate decay rate for the velocity if we choose equation (3.21). There is no 
correction on the density but the Fourier mode bo is 0 from assumption 3.1. 

We can see that if we use (3.22), then this choice ensures non-negativity of all 
Fourier coefficients of the convolution kernels, and from (3.15), we see that —Xdk is 
at least d. For small modes, the diffusion process is not large enough to ensure the 
decay rate, so we need to add the feedback term. For larger modes, diffusion will be 
enough, and there is no need to add the feedback (but one can still add a feedback 
term, and the decay rate will become larger than the specified rate for such modes). 
Note that even if we drop the max{0;.} in (3.22), it still ensures that the decay rate 
will be at least d, but niay become negative for large modes, and this can lead to 
numerical instabilities. 


Then, recalling equation (3.18), we set (ppk as in (3.23). This choice of (ppk 
ensures that < 0 , and the decay rate corresponding to eigenvalues X±k is then 
(Cl -I- C2)|fcp -I- (fiuk 


exactly 


which is always larger than (or equal to) d, because of 


the choice of ipuk in equation (3.22). 

The max{0,.} in equation (3.23) is also optional, since even without the max, the 
discriminant will be negative, but the non-negativity of tppk will avoid numerical 
instabilities. 


For large modes, namely |A:| > — and |A:| > ./ ^ , ip^k = 0, and then there 

is no issue in considering the inverse Fourier transform and defining a convolution 
kernel <pu{x) with these coefficients. 

With the above choices, only finitely many Fourier modes of <pu{x) and <Pp{x) are 
nonzero and hence the convolution kernels exist, which proves the result. □ 


Note that we could use the same definitions from equations (3.22)-(3.23) for larger 
modes (|A:| > K), but then the Fourier series for p>p{x) does not converge, as (ppk = 
Odfcp). So we need to truncate the series, and set (ppk = 0 for |A:| > K. 

We also note that for large modes for which both (ppk and (puk are set to 0, the 
observer is simply a solution of the equation without any forcing term, and the largest 
eigenvalue is A+fe (both Xdk and X-k —>■ —oo when |A:| —>■ -l-oo). From (3.15), the decay 

7 

rate is then equal to ——— = asymptotically for \k\ -|-oo. 

Cl + C 2 A “h 2fi 


3.5. Observers with density observations. We have only considered obser¬ 
vations of the velocity field so far in this paper. In this subsection, we will construct 
the observers for the case when the density is observed instead of velocity. 

Since observations of the scalar field (density) naturally gives less information 
than those of a vector field (velocity) in any dimension greater than one, we will see 
that, as expected, we will not be able to control all the modes in the observer. In 
particular, we will be unable to modify the convergence rates for the eigenvectors 
corresponding to the diffusion eigenvalues but we will be able to get arbitrarily high 
convergence rates for the eigenvectors corresponding to eigenvalues A±. Since the 
details of the calculations are very similar to those presented above, we will only 
present a sketch of the results. 

When we have density observations p{t, x) instead of velocity u{t, x), the observer 
is exactly as in the observer equations (2.8) with u being replaced by p. Then the 
equations for errors r = p — p and v = u — u will be as follows. These are the same 
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as equations (3.1) with r replacing v in the feedback terms. 

(3.25) r, + p„V.. = -^,.i),r, 

Povt = pAv + (A + ^)V(V • v) - 7 Pq Vr -tpu* DuV . 

We see that because of the symmetries of the equations, the natural choice of the 
differential operators Dp and Du with minimal derivatives is 

(3.26) Dp{f) = f and .D„(/) = qpj"V/. 


Assumption 3.7. We suppose that the mean value of v{t,x) is equal to 0, i.e. 
the mean values of the initial conditions of u and u are the same. 

Again, this is a strong assumption similar to assumption 3.1: contrary to the 
density mean, the velocity mean cannot be corrected from (3.25)-(3.26). 

Remark 3.8. As in the case of velocity observations, we can eliminate v in order 
to obtain an equation for the density alone. This can help understanding the choice 
of the differential operators in the feedback terms. From the first equation of (3.25), 
we can extract pfS/ ■ v = —rt — ifp * Dpr and then we need to take the divergence of 
the second equation: 

PqV • ut = /iV • Au + (A + p)V • (V(V • v)) - • (Vr) - ■ {Dur). 

By choosing Dp and Du as in (3.26), we get the damped wave equation for r{x,t) 


(3.27) p^rtt = (A + 2p)Art + jPoAr + (A + 2p)'ipp * Ar + jpfitpu * Ar - pQtjjp * rt , 

which again shows that the choices of Du and Dp are quite natural to control the terms 
in the above equation. 

We can then see that if we are in one dimension n = 1, equations (3.4) and (3.27) 
are identical with ipu = ifp, and (pp = ■!/)u + ((2/i+A)'0p)/(7p2). Thus in one dimension, 
observations of the density and those of velocity give the same results for rates of 
convergence of the observer solution. 

Introducing the Fourier decomposition for (r, v) exactly like in equations (3.6) 
and the vector yk of the Fourier coefficients, we will get a system 2 /(,(t) = Mkyk{t), 
where Mk is the following (n + 1) x (n + 1) complex matrix: 


(3.28) 


Mk = 


-ippk Cik"^ 

-C3(l + ^„fe)fc -Ci|fcp/„ - C2fc I 


po*27r, C5 = C3 X C4 = 'yPq ^47r^. 

We also recall that {/c^,l<i<n—l}isa basis of the set of vectors orthogonal to 
k in M" if fc 7 ^ 0. We have then the following result: 

Proposition 3.9. 

• For k ^ 0, the eigenvalues of matrix Mk are: 

(3.30) 

Adfc = -ci|fcp, 

- ((Cl + C 2 )|fcp + l/jpk) ± ((ci + C2)|fc|2 - tjjpkf - 4 c5(1 + V'«fc)|fcp 
A±fe =- 


We recall the values of the constants: 
(3.29) 

/i47r^ (A + /i)47r^ 


Cl = 


Po 


C2 = 


Po 


; C3 = 7Po ^*27r, C4 = - 


2 
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The eigenvalue Xdk ho-s a multiplicity n—1, with eigenvectors 

eigenvalues X±k both have a multiplicity 1, with eigenvectors ^ 
For k = 0, the eigenvalues of matrix Mq are: 


The 


0 

C4|fcp 

(A±fc +i>pk)k 


(3.31) 


Xdo = 0, A_o = —4’po- 


Xdo now has multiplicity n, with associated eigenvectors: 
i < n} beim 
eigenvector 


, {ei,l < 


i < n} being a basis o/R", and A_o has multiplicity 1, with the associated 

1 


The proof of this result is straightforward, one simply has to multiply these eigen¬ 
vectors by Mfe. 

Note that for k = 0, Xdo cannot be controlled, but the associated eigenvectors 
correspond to constant velocity solutions. So from assumption 3.7, this case does not 
appear. But in this case A_o can be controlled, and it corresponds to constant density 
solutions. 

For fc ^ 0, the discriminant is 


(3.32) Afc = ((ci -I- C 2 )\k\'^ - ippk)'^ - 4 c5(1 -I- '0ufc)|fc|^ 


and we still use the notation \/= iyj—ls.k if < 0. 

Note that Xdk is independent of ijjp and ipu and thus cannot be controlled by any 
choice of tpp or Exactly in parallel with the previous case (velocity observations) 
and the choice presented in equations (3.22)-(3.23), we have the following result: 

Theorem 3.10. We suppose that assumption 3.7 holds. For any d > 0, for any 
K > 0, one can find ipp{x) and ipuix) such that the maximal decay rate of the errors 
r{t,x) and v{t,x), solutions of (3.25), towards 0 is at least min{d,cijfcp} for any 
Fourier mode k such that |A:| < K. The following values can be chosen, with ci, C 2 
and C 5 given by (3.29); 


(3.33) tjjpo = d, tpuo = 0 , 

(3.34) = max{0; 2d — (ci-I-C2)|fc|^)} , 0<|fc|<iF, 

I j ^ {{ci + C2)\k\^ - i^pkY \ 

(3.35) V'«fc=max<0;--^f’ 0<|fc|<Ar, 

(3.36) ^/;„fe = 0, iljpk = 0, \k\>K. 


The max(., 0) in (3.35) is there to ensure positiveness of all Fourier coefficients of ipu- 
It could be relaxed as it is only useful for small modes (the other term grows as |fcp). 
Proof. The choice of ipuk ensures that A^ < 0 and the real part of X±k is then 
({ci + C2)\k\'^ + fjpk) 

- - - - -. Then, using (3.34), the decay rate is then at least d for all k 

with \k\ < K. n 

We note that the modes that we cannot control to have a pre-specified decay 
rate are the incompressible ones since the eigenvector 


kf 


satisfies k ■ Vk = 0, 


corresponding to V • u = 0. This is expected since the density observations cannot 
give any information about the “constant density” incompressible flow. 
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3.6. Remarks and comparison with the nudging feedback. We now com¬ 
pare the result presented in the previous sections with what can be done with the 
standard nudging observer, for which only the observed variable is corrected. We will 
use a = (ci -I- C 2 )|fcp, b = 2,y^\k\ (see (3.29)) in the discussion in the rest of this 
section. 

Nudging with velocity observations. In this case, there is no feedback on the den¬ 
sity: ifp = 0, and the differential operator on the velocity Du is simply the identity 
(or more physically Du{f) = pof as in (3.3)). Then equations (3.1) become 


(3.37) 


rt -I- poV • u = 0 , 

povt = pAv -I- (A -I- /r)V(V • v) - ^pl~^S/r - pQipu * v . 


Let us fix a: = ipuk for this paragraph. We can write the eigenvalues A± from equa¬ 
tion (3.15) as 


A± 


(a -I- a;) ± \/(a -|- a;)^ — 6 ^ 


From this expression, we can find optimal decay rates as follows. 

• If a < 6 , i.e, for |fc| < —^± is 

maximized for x = b — a, and for this choice the decay rate is 6/2 = y^|A:| 
which is greater than the decay rate without nudging (x = 0) but it cannot 
be made arbitrarily large as in the full state observer (see theorem 3.6). 

• If a > 6 , i.e., for |fc| > again see that the real part of —X± is 

maximized for x = 0 and this simply gives the decay rate without nudging. 
Again, the full state observer presented above is better then simple nudging 
since we can obtain arbitrarily large decay rates up to |fc| < K for any fixed 
K. 

Nudging with density observations. In this case, there is no feedback term on the 
velocity: ipu = 0, and the differential operator on the density Dp is still the identity, 
as in (3.26). The system (3.25) now becomes 


(3.38) = 

poVt = pAv -I- (A -I- /r)V(V • v) - ypp Vr . 

Let us fix a: = il^pk for this paragraph. We can write the eigenvalues X± from equa¬ 
tion (3.30) as 


A± 


(a -I- a;) ± \/{a — xY — 6 ^ 


From this expression, we can find optimal decay rates by maximizing the real part of 
the above expression as a function of x in exactly the same manner as the discussion 
in the previous paragraph. In this case, it is easy to see that the real part of —A± 
is maximized for x = b -\- a, and for this choice the decay rate is min{a -I- 6/2,ci|fcp} 
which is greater than the decay rate without nudging (x = 0 ) but the first factor in 
the min cannot be made arbitrarily large, and the decay rate may be smaller than 
the one for the full state observer (see theorem 3.10). 

Thus we see that in both cases of either velocity or density observations, the full 
state observers previously presented perform better than simple nudging, even with 
an optimal choice of the nudging convolution kernel. 
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4. Additional study of the observer in dimension 1. In this section, we 
study the same observer as in the above discussion, but in one dimension. There are 
a few simplification which we point out below. We extend the study of the above 
observer to two special cases in one dimension. In subsection 4.2, we discuss the 
case when the equation contains an unknown forcing term while in subsection 4.3, we 
consider the case when the observations are over a subdomain [0,L] with L < 1. 

4.1. Remarks about decay rates in one dimension. Essentially all the 
results in the previous section are also applicable in one dimension with the only 
difference being the absence of certain eigenvalues, as we discuss in this subsection. 
Noting that fj,Au + (A + /i)V(V • u) = in one dimension with v = 2fi + X, we can 

simplify equation (3.1) to the following: 

(4.1) rt + poVx = -po^fp * Vx , Povt = vvxx - ip1~^r^ - popu * v , 
whereas equation (3.4) simplifies to 

(4.2) povtt = vvtxx + IpIv^x + 'yp^^Pp * Vxx - Pq^Pu * Vt ■ 

From this, we can see that the matrices and iV^ from equations (3.10) and (3.13), 
respectively, will remain the same, noting that k ® k = |fcp = k"^. These 2x2 
matrices will not have eigenvalues Aofc or Xdk but only two eigenvalues X±k as given 
in equation (3.15) with the same eigenvectors as given before. 

Hence the main result in theorem 3.6 also applies in one dimension, with the 
following modification to the choice of (puk- 

(Puk = max{0; 2d — (ci + C 2 )k^)} = max <0;2d - 

I Po 

and the same choice of Ppk as in equation (3.23): 



(fpk = max • 


V ((ci 


C2) 


4c5|fc|" 


— 1 > = max < 0; 


, PO 


^uk 


167r^7PQ ^k'^ 


- 1 


Further, in the case of density observations, the result in theorem 3.10 is even stronger, 
in the sense that with the choice of il)pk and '(/'ufe of equations (3.34)-(3.35), the decay 
rate is at least d, since there is no non-trivial incompressible flow field in one dimension. 

4.2. Unknown forcing term. We now consider observers for the linearized one 
dimensional equations with a forcing term in the velocity equation. For simplicity, we 
also assume that po = 1, and as earlier, without loss of generality, we assume uq = 0. 

(4.3) p^+Ux = 0, Ut = i'Uxx-lPx- f{x,t). 

If this forcing term f(t, x) is known, then we also add it to the velocity equation of the 
observer (2.8). Then, by considering the difference between the observer and original 
equations, the forcing term disappears and we are still considering equation (4.1) for 
the error. 

So we now assume that the forcing term is unknown. In this case, we cannot add 
it inside the observer equation. So the observer equation remains unchanged, and 







15 


then, the difference between the reference velocity and the observer velocity satisfies 
equation (4.2) with a forcing term: 

(4.4) Vtt = VVtxx + IVxx + * Vxx - ‘^u*Vt + ft- 

We now adapt the spectral analysis. We assume that the mean of / is 0 (no bias 
in the forcing), and then the Fourier decomposition of / is: 

k^O 

Then from (3.12), we obtain the following new equation for the mode of the error: 

(4.5) + {ipuk + 4i/fc^7r^) a'^(t) + 4fcV^7(l + (ppk) ak{t) = 

We just need to find a particular solution to this equation, and add it to the general 
solution that we found in section 3.3. 

We consider a very simple case (for clarity reasons and for a better understanding 
of the phenomenon), where the time dependence of the forcing is a sine (or cosine) 
function: 


(4.6) 


fk{t) = CfcSin(2a;fc7rt), 


where Wk is the frequency of the forcing oscillation of mode k. Defining ak = ^uk + 
Ai'k^TT^, and Pk = 4fc^7r^7(l + (ppk), equation (4.5) becomes: 

(4.7) alit) -j-aka'kit) PkO-kit) = 2cfeWfe7r cos(2wfc7rt). 

A particular solution is then given by 

(4.8) ak{t) = Ak cos(2a;fc7rt) + Bk sin(2a;fc7rt). 


The constants Ak and Bk are solution of the following linear system: 

Afc(/3fc - 4u;fc7r2) + Bk{2uJkTrak) = 2cfcWfc7r, 

Afc(-2wfc7rQ:fc) + BkiPk - 4Wfc7r2) = 0. 


Then, we get: 

(4.9) 

(4.10) 


Ak — 2ckUJk'^ 


_ Pk - 4cJfc7r^_ 

{Pk - + (2a;fc7rafc)2 ’ 


Bk 


^ _ 2u)k'^kXk 

{Pk-Aujl'K^f + {2u:k^akY' 


The amplitude of the particular solution is then given by 


(4.11) 



2CkUJk'^ 

-YiPk - + {2uJkTTakY 


Increasing ippk or ipuk (or both) will make Pk or ak (or both) increase, and then 
Dk will decrease. This means that we can make Dk become as small as we want 
and the observer will converge towards the true state. But of course, the numerical 
performance of the observer is severely degraded in comparison with previous cases. 
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as it is usually not possible to consider extremely high values of feedback coefficients 
from a numerical point of view. 

Concerning the density, adapting the previous calculations on the decrease of p 
(knowing the decrease of u), we get the following amplitude of the particular solution: 

(4.12) Ek = —(l + ippk)-Dk, 

which means that one can adapt the limit amplitude of p by changing the values of 
ippk- Theoretically, choosing ippk = — 1 has the effect of completely removing the 
influence of the forcing term on the density, but of course, it is not numerically stable 
or physically consistent to consider negative feedback coefficients. 

4.3. Heuristic in the case of observations over a part of the domain. 

In the case of observations over only a part of the domain, we define the observation 
domain O = [0, L] with L < 1. In this case, in order to proceed with Fourier analysis, 
we will need to find the Fourier transform of r;lQ because the feedback terms ip * v 
will be replaced with (p * (ulo). The Fourier series for 1 q is 

• r, 7 1 ~ cos 2TrkL 

sin ZTTKX H- cos ZTTKX . 

k 

Thus we see that the Fourier series of ulo will have Fourier components for all k even 
in the case when v has just a single Fourier mode. 

Such coupling of Fourier modes will numerically slightly degrade the performance 
of the observer, as at any time, some energy will be transferred between different 
Fourier modes, as in the full nonlinear model. 

If we rewrite equation (4.2) (for simplicity with po = 1), we get: 

(4.14) vtt = vvtxx +'jv XX * {,'^XX In) — (vtin) ■ 


(4.13) 


In — T + 


fc=i 


sin27rfcL 


Then, for simplicity reasons, assuming v only has one single Fourier mode k, equa¬ 
tion (4.5) rewrites: 

(4.15) a'l{t) + {Lipuk + 4z/fc^7r^) a'^{t) + 4A:^7r^7(l -|- Lppk) ak{t) = 0. 

Indeed, only the 0th order term {L, from (4.13)) in the Fourier decomposition of 1 q 
will be kept through the convolution with v (or one of its derivatives). So, the decay 
rate becomes 


(4.16) 


<fukL + ivk^TT^ 
2 


Of course, this is an approximation, as even if v only has a single Fourier mode at 
time t = 0, the convolution with a characteristic function leads to mode mixing, as in 
a nonlinear situation. But we assume here that most of the energy is along mode k (if 
only this mode is present at the initial time), which will be confirmed by numerical 
experiments in next section. 


5. Numerical experiments on the ID compressible Navier-Stokes ob¬ 
server. In this section, we report some of the numerical investigations in one dimen¬ 
sion that illustrate the linear theory we discussed above. We also present results of 
using the same observer as in the linear case for two additional scenarios, namely, (i) 
the observer when the velocity is observed only over a subinterval of the domain in 
section 5.4, and (ii) in section 5.5, observer for the fully nonlinear equations. 
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Fig. 1. The norm of the difference between the observer {p,u) and the solution {p,u) versus 
time. Solid and dotted lines are the errors in p and u, respectively. The left panel is for fixed tpp = 0 
with varying tpu while the right panel is for fixed ipu = 20 with varying (pp. 


5.1. Numerical configuration. The space domain is [0, 1] and we still assume 
periodic boundary conditions. The discretization involves 10^ grid points, with a 
step Ax = 10“^. We consider a time step At = 10“^. We also experimented with 
increasing the spatial resolution (and correspondingly decreasing the time step), but 
the results are almost identical and not presented here. The adiabatic exponent is 
set to 7 = 1.4, and the diffusion is set to ly = 5.10“^ (except in section 5.5). The 
numerical code uses a conservation form of the compressible Navier-Stokes system, 
with p and pu as variables. A hnite volume scheme is used, in which the inviscid 
flux is computed using an approximate Riemann solver (e.g. VFRoe scheme). Time 
integration scheme is a third order explicit Runge-Kutta scheme, where the time step 
is chosen based on a CFL condition. 

In order to reproduce a quasi-linear situation, we consider the true (observed) 
solution p(t,x) = 1 and u{t,x) = 0 while the initial conditions for the observer are 
set to: 


(5.1) 


piix) = 1 + 5.10 ^ sin(27rfcx), u/(x) = 5.10 ^ sin(27rA:x), 


so that the mean values of p and u are pi = \ and uj = 0 respectively, and where k 
is a given mode, usually the first one (fc = 1, unless differently specified). 

We first look at the solution without any feedback term, see equation (2.8) with 
<fu = <fp = 0- Figure 1 (solid and dashed curves, left panel) shows the evolution (in 
log scale) of the norm of the difference between the observer (p, u) and the solution 
(p, u). As there is no feedback, all the Fourier coefficients cppk and ipuk are equal to 0, 
and then from equation (3.18), the discriminant is A ~ —217.2, and the theoretical 
decay rate (only due to diffusion) is given by equation (3.15): dth = 0.987. Also the 

Att 

oscillation period can be computed from (3.15) and (3.18): ojth = —j= — 0.85. 

Numerically, the slope of the solid curve in the left panel of figure 1 gives a 
numerical decay rate c?num = 0.980. Note that the figures show the errors to base 
10, hence the slope of the semilog plot is dnum/log(lO) = 0.426. The numerical 
oscillation period is approximately uinum = 0.852. Note that one period corresponds 
to two oscillations on the figure for the norm of the cosine. This excellent agreement 
between theoretical and numerical values can be reproduced for other modes and 
other values of the parameters. 
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Theoretical 
decay rate 

Numerical 
decay rate 

Theoretical 
oscillation period 

Numerical 
oscillation period 

0 

0.987 

0.980(0.426) 

0.85 

0.86 

0.1 

1.037 

1.032 

0.85 

0.85 

0.5 

1.237 

1.237 

0.86 

0.85 

1 

1.487 

1.486 

0.86 

0.87 

5 

3.487 

3.485 

0.97 

0.98 

10 

5.987 

6.012(2.61) 

1.42 

1.38 

12.895 

7.434 

6.590 

81.0 

— 

15 

4.393 

4.364(1.895) 

- 

- 

20 

2.897 

2.861 

- 

- 


Table 1 


Theoretical and numerical decay rates and oscillation periods for several values of ipu (k = 1, 
ifp = Oj. The values in parenthesis give the numerical decay rates in base-10, in order to compare 
with slopes of lines in figure 1 


5.2. Simple nudging observer. We now consider the nudging framework (see 
section 3.6 and equation (3.37)), and we hence suppose that there is some feedback 
only in the velocity equation. We then first let (fp = Q, and only modify the values 
of ipu- This simulates the nudging, or asymptotic observer: as only the velocity is 
measured, only the velocity is corrected in the observer system. Table 1 shows the 
theoretical and numerical decay rates and oscillation periods for several values of ifiui- 

The first remark is that the numerical results perfectly match the theoretical 
results, except for the particular value of (pui = 12.895. In this case, the numerical 
decay rate is slightly smaller than the theoretical one. Also, no oscillations can be 
seen on the results, which is reasonably in agreement with a theoretical period of 81 
which will be impossible to see with a final time of T = 5. 

As ipui increases, the decay rate increases, until ifiui reaches dTr^/y—4z/7r^ ~ 12.895 
(see remark in section 3.6), for which the discriminant is equal to 0, and then positive 
for increasing ipui- The corresponding optimal decay rate is dth = — 7.434 (see 

section 3.6). We can see that the decay rate then decreases, as the discriminant takes 
larger positive values, so that one of the two eigenvalues gets closer to 0. 

One can see on figure 1 (left panel, dotted and dash-dotted curves) that the 
error decreases much stronger than the case of no feedback (figure 1, left panel, solid 
curve). We can also see that with increasing kpu, the period of oscillations increases 
and eventually there are no oscillations (discriminant is positive). It confirms that 
the decay rate decreases if fiux is increased too much. 

Similar results have been observed for other values of the diffusion v, and for other 
modes fc, and in each case, these numerical results match well with the theoretical 
predictions. 

5.3. Results on the full state observer. We now use the full state observer, 
with an additional feedback term in the density equation (see equation (3.5)). We 
refer here to theoretical results from sections 3.3 and 3.4. We first set (fui = 20, for 
which the discriminant is positive, the largest eigenvalue gets closer to 0, and then 
the decay rate becomes non optimal. 

Table 2 shows the theoretical and numerical decay rates (and oscillation periods) 
for several values of (ppi. As the discriminant is positive when (fui = 20 and (ppi = 0, 
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‘Ppl 

Theoretical 
decay rate 

Numerical 
decay rate 

Theoretical 
oscillation period 

Numerical 
oscillation period 

0 

2.897 

2.861 

— 

— 

0.5 

4.838 

4.790(2.080) 

- 

- 

1.184 

10.94 

9.74 



1.5 

10.98 

11.03(4.791) 

1.50 

1.57 

5 

10.99 

11.03 

0.43 

0.43 

10 

10.99 

11.06(4.802) 

0.28 

0.28 


Table 2 


Theoretical and numerical decay rates and oscillation periods for several values of y>p with fixed 
if>u = 20 for the k = 1 mode. (The values in parenthesis are again decay rates in base-10 for 
comparison with figure 1.) 


equation (3.23) gives the theoretical value for which the discriminant comes back to 

+4z^7r^)^ 

negative values: ippi = --—^-1 ~ 1.184. 

16771^ 

We observe numerically that the decay rate is very similar to the theoretical rate, 
and we clearly observe the transition from positive to negative discriminant with the 
stabilization of the decay rate, and the apparition of oscillations. Increasing ippi to 
a much larger value than the optimum given by (3.23) is not necessary, as the decay 
rate does not increase, and the period of oscillations increases quite quickly. This is 
clearly seen from the plots in the right panel of figure 1. 

As we have seen, by adding the derivative of the velocity as a feedback to the 
density equation, we were able to significantly increase the decay rate of the error, in 
comparison with only a feedback in the velocity equation. 


5.4. Observers with observations over a part of the domain. In this 
section we present the numerical results of the full state observer, see equation (2.5), 
but with observations over only a part of the domain. The left panels of figure 2 shows 
the decay rate of the Lp' error between the observer (p, it) and the actual solution (p, u). 

We see that, as expected from section 4.3, the rate of decay is smaller for smaller 
observational intervals. We also see that the error in velocity decreases linearly (with 
oscillations) whereas the error in density saturates at a fairly high but constant value. 
This is because the solution of the observer equation converges to a solution with 
ti = 0 but with p = p„ where Pn ^ Pi - i-e. the observer density is shifted by an 
amount which compensates for the initial discrepancy between the mean of p and the 
mean of p. This is not surprising since the original equations themselves are invariant 
under the constant shift in density. 

In order to overcome this problem in the case of observations over partial domain, 
we propose the following modification of the feedback terms in equation (2.5): 


(5.2) T’p(p, It, u) = ipp* Dp [la(M -u)- {u- it)], 

where (/) indicates average of / over the interval [0, L]. This ensures that the average 
of the feedback term is zero and hence the equilibrium solution of this equation also 
has mean zero. Note that in the case of L = 1, i.e., the case of full observations, 
this average is just zero and the feedback in equation (5.2) is identical to that in 
equation (2.5). The errors obtained by using this new observer are shown in the 
right upper panel of figure 2. We clearly see that the observer p now approaches the 
solution p and the error decreases as expected. 
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Fig. 2. norm of the difference between the observer p,u and the solution p,u, in the case 
of observations over the sub-interval [0, L] for L = 0.1, 0.3, 0.5, 0.8,1, all of them with strength of 
feedback terms being tpp = 0.5 and ifu = 10. Left panels are for feedback from equation (2.5) while 
the top right panel is for feedback from equation (5.2). The error in velocity for the feedback from 
equation (5.2) is almost identical to the bottom left panel and hence not shown. The bottom right 
panel shows the observer solutions at time t = 0.16. Note that “observer 1” refers to (2.5) while 
“observer 2” refers to (5.2). 


In order to clearly see the effect of the observer, the lower right panel of figure 2 
shows the actual observer solutions for the case when If = [0,0.3]. They clearly 
show the effect of incorporating the observations, and also the difference between the 
observer with the incorrect mean and the one with correct mean. The effect is of 
course more pronounced on the density than on the velocity. 

In this case the discriminant of the decay rate is clearly negative (as evidenced by 
oscillations in the plot for errors) but we can calculate the decay rate. Figure 3 shows 
the decay rate of the observer as a function of the length of the interval over which 
the velocity is observed, for the case when (pp = 0.5, (pu = 10 for the k = 1 mode. 
We notice that this is pretty close to a straight line and the best fit line is given as 
follows: 


(5.3) 


D = 4.99L + 0.24 . 


Comparing this with the rate given in equation (4.16), we see that the last equality is 
a reasonable approximation. Thus even though we cannot calculate the exact decay 
rate (see section 4.3), a reasonable estimate can be obtained by using formula (4.16): 


(5.4) 


D 


PuL + 4z/7r^fc^ 
2 


= 5T +0.987 


in experiments corresponding to figure 3 ((/?„ = 10). 





















































21 



0.2 0.4 0.6 0.8 1 

Width of observed interval (L) 


=0.5,=10.0,fc = l 


•—• Density error decay rate 
Velocity error decay rate 

— Best fit line (density) 

— Best fit line (velocity) 

. L +4i^^ )/2 


=0.5,=10.0,fc = l 


•—• Density error decay rate 
Velocity error decay rate 
- — Best fit line (density) 

— Best fit line (velocity) 

. -L +4i/7r^ )/2 
























0.2 0.4 0.6 0.8 

Width of observed interval (L) 


Fig. 3. The decay rates for the velocity and the density, as a function of the length of the 
observation interval (quasi-linear case in the left panel and fully nonlinear case in the right panel). 
Note that the slope of the best fit line ^ 4.5 whereas the slope of best fit line by taking only the points 
for L = 0.1, 0.2, 0.3, 0.4, 0.9,1.0 is ^ 5.7. 

This figure shows that the global slope of the numerical decay rates versus the 
length of the observation interval is close to ^ = 5. There is also a good agree¬ 
ment between numerical decay rates and the approximated theoretical ones (see equa¬ 
tion (4.16)) when L is either small (almost no observations) or large (almost all the 
domain is observed), while the numerical decay rates are degraded when only half 
the domain is observed. When L Ri 0.5, previous spectral studies show that mixing 
of Fourier modes has a higher effect than when L « 0 or 1. Convolution with the 
characteristic function leads to non negligible transfers of energy from mode k to other 
modes, and then to a smaller decay rate. 

But figure 3 also shows that even with a small observed subdomain, our observer 
is still very efficient and we can still control the decay rate by increasing . We note 
that we have not considered the problem of placement of the observations, which is 
certainly important and has a huge influence on the decay rate [29, 19]. 

5.5. Observer for nonlinear Navier-Stokes system. In this section, we will 
report the numerical results of using the full state observer for the nonlinear system of 
equations, i.e. from equations (2.3), with the feedback terms as in equations (2.5), in 
one dimensional system. Note that we do not have theoretical estimates of the decay 
rates towards the equilibrium solution. The aim of this section is to understand the 
efficiency of the above nonlinear observer. 

Firstly, we consider observations of the equilibrium solution, so that we are essen¬ 
tially studying the decay of the observer towards the equilibrium. The left panel of fig¬ 
ure 4 shows the decay rates of the observer solution with the following two initial con¬ 
ditions, with the feedback term set to zero, and for the case with {ipp, tpu) = (0.2,10): 

Pi{x) =1-1- 5.10“^ sin(27rA:a;), ui{x) = 5.10“^ sin(27rA:a;), 

Note that the perturbation strength is 10 times larger than the “quasi-linear” case 
with the perturbation in equation (5.1). We consider two cases, k = 1 and k = 3. 
The cases with other values of fe > 1 show a very similar behavior to the case fc = 3 
and hence is not shown here explicitly. For reference, we have also plotted the decay 
rates of the linear case as well. 

The right panel of that figure shows the actual observer solutions for the k = 1 
case. Note the different scales for the nonlinear (left axis) and linear (right axis) 
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Fig. 4. The decay rates of the nonlinear observer compared with that of the linear observer 
(top left). The actual observer solutions are also shown (k = 1 in the top right, k = 3 in the bottom 
right). The decay rates for the nonlinear observer with observations over partial domain is in the 
bottom left. 


regimes. We clearly see that with perturbation strength of 5.10“^, even though the 
initial condition only has k = 1 mode present, at time t = 0.1, higher modes are excited 
(dotted line, with no feedback). When the observer feedback is added, the solution 
very quickly decays before higher modes are excited, bringing the perturbation to 
the level where the linear theory is a good approximation and thus the decay rate is 
identical to that of the linear case of perturbation strength 5.10“^. 

The case for fc = 3 (and indeed all higher modes with A: > 1) is quite interesting. 
In this case, the initial condition of fc = 3 excites the fc = 1 mode. Thus even when 
the observer feedback is added, the decay rate is not the same as the linear fc = 3 
decay rate but rather it is closer to the linear fc = 1 rate. The same behavior is seen 
for other modes with fc > 1. Figure 4 shows example of this decay for fc = 3 and 
fc = 6, along with the actual observers. It is quite clear that in the quasi-linear case, 
modes other than the one contained in the initial condition are not excited while in 
the nonlinear case, they are quite clearly excited. 

Finally, we also performed numerical experiments, still with the full state observer 
on the nonlinear system, but now with partial observations over various domain sizes, 
as discussed in section 4.3. Quite surprisingly, the behavior in the nonlinear case is 
very similar to the linear case: the decay rate is close to being linear in the size of the 
domain - see the right panel of figure 3 and the bottom left panel of figure 4. 

6. Conclusion. In this paper, we were interested in observer design for a viscous 
compressible Navier-Stokes equations. Thanks to intrinsic properties of the system 
(symmetries), we were able to design observers in order to reconstruct the full solution 
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(both velocity and density) when only one variable (either velocity or density) is 
observed, even partially. 

A spectral study showed that for a tangent linear system (linearized around the 
equilibrium state of constant velocity and density), we can prove the convergence of 
the observer towards the solution, and we can control the decay rate of the error, with 
explicit formulas for the feedback coefficients as functions of the desired decay rate. 

Numerical experiments in one dimension are in perfect agreement with theory 
in the linearized situation: we can obtain any decay rate by increasing the observer 
coefficients. Numerical experiments also show that our observer is still very efficient 
in the case of observations over partial domain, and also in the full nonlinear case. 

The application of this kind of observer in other specific cases such as full prim¬ 
itive equations of the ocean or the atmosphere, will be a natural extension of this 
work. Observers for fluid equations coupled to other quantities such as temperature 
or salinity with observations of these fields instead of velocity observations will also be 
an interesting extension that will be of great interest in practice since these types of 
measurements are more common than measurements of the velocity field. It will also 
be quite challenging and interesting for practical applications to consider observers in 
the case when observations are available discretely in time and/or in space. 
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